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Oooan-lhve  Record  Analysis  - Ordinate  Distribution  tad  Way  Heights 


ABSTRACT 


She  methods  of  Analysing  ooesn  vara  reoorda  into  ordinates  and  into  waves 
are  oompared*  It  is  suggested  that  aui  table  wars  spectra  will  permit  many 
of  the  numerical  res  alts  of  one  procedure  to  be  obtained  from  the  other* 

The  Rayleigh  distribution  is  introduced  as  that  of  the  envelope  of  a 
Gaussian  randan  process,  whioh  approximates  the  peaks  and  trouts  of  the 
prooess  ourve  in  the  oase  of  narrow-band  speotra.  This  distribution  is 
shorn  to  resemble  ware -height  distributions  olosely  in  regard  to  shape, 
dependence  on  a single  parameter,  and  values  of  various  parameter  ratios* 

The  single  parameter  may  be  interpreted  as  the  dispersion  (r*m.s.  deviation) 
of  the  ordinate  distribution,  or  the  total  power  in  the  wave  spectrum* 

The  ordinates  of  twenty-one  time  histories  analysed  on  the  Wave -Ordinate 
Distribution  Analyser  su*e  found  to  follow  the  Gaussian  distribution  olosely 
enough  to  permit  ready  graphioal  estimation  of  the  ordinate  dispersion* 

The  mean  wave  heigit  is  found  to  be  olosely  proportional  to  the  .latter 
quantity,  the  proportionality  fhotor  being  about  10jC  1ms  than  that  for 
the  mean  envelope  height*  The  so&cter  around  the  straight  line  of  pro- 
portionality is  small  enough  to  indicate  that  measurement  of  mean  wave 
height  may  for  many  purposes  be  replaoed  by  measur ament  of  wave-ordinate 
dispersion,  whioh  is  more  suited  to  analog  computation* 

The  present  data  are  used  to  illustrate  simplified  wave-ordinate  distri- 
bution analyses  yielding  reasonably  ao curate  estimates  of  ordinate  dispersion 
and,  hence,  mean  wavs  heigit*  There  is  suggested  an  ocean-wave  measurement 
program  based  on  these  results  designed  to  provide  a neerly-oontinuous  flow 
of  immediate  ware -height  information  over  a recording  period  extending  over 
several  seasons* 


le  Introduction 


Visual  observation  of  the  surface  of  the  rea  has  in  many  — ses  led  to 
its  analysis  into  a series  of  orests  and  troughs.  Jbr  a mildly  irregular 
s—  r— orded  at  a point,  such  features  appear  on  a time  history  of  the 
surface  or  of  the  subaurfaoe  pressure  as  well-defined  prominences  and  de- 
pressions on  a continuous  ourve,  The  more  regular  the  fluctuation,  the 
more  nearly  is  this  type  of  decoription  adequate,  and  in  the  limiting  — se 
of  perfect  regularity,  knowledge  of  the  height  of  any  one  wave  and  the  period 
between  any  two  euooeesive  wave  or— ts  ie  euffloient  to  oharaoterise  the  re- 
sulting sinusoidal  wave*  The  more  irregular  the  time  history,  on  the  other 
hand,  the  less  complete  will  be  the  information  contained  in  the  results  of 
an  enumeration  of  "wave  heights*  and  “wavs  periods"*  In  faot,  it  has  been 
found  that  for  typioal  records  of  — earn  swell  the  entire  statistical  distri- 
bution of  estimated  surftoe  wave  heights  over  a twenty-minute  interval 
yields  essentially  only  one  number  for  characterising  the  wave  record* 

In  view  of  the  discussion  below,  this  may  be  taken  as  evidenoe  that  such 
wave  systems  have  speotra  with  frequency  bandwldths  whioh  are  narrow 
relative  to  the  oentral  frequenoy  of  the  band.  It  is  reasonable,  of 
course,  to  —ex me  that  the  speotra  of  wave  reoordj  having  oiailar  geo- 

have  important  oamacn  features  characterising  their  shape* 


2. 


8uoh  theories  suggest  that  oertain  statistical  information  about  ware 
heights  and  periods  nay  eventually  contribute  substantially  to  sub  ade- 
quate description  of  the  state  of  the  sea* 

However,  until  more  is  learned  about  the  form  of  typical  ocean-wave  spectra, 
it  is  natural  to  direot  attention  to  the  continuous  aspeots  of  the  wa.ve- 
reoord  ourre,  rather  than  to  its  discrete  wave-by-wave  oharaoteristios* 

The  oonoept  of  a Iburier  speotrum'®'*  then  beoames  available,  together  with 
the  results  cf  other  linear  operations  performable  on  the  ordinates  of  the 
time-history  funotion,  scene  of  whioh  may  be  conveniently  arrived  at  by 
analog  computation. 


i 

i 
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The  purpose  of  the  present  report  is  to  present  certain  experimental  results 
which  relate  the  ordinate-wise  and  wave-wise  approaohes  to  Idle  analysis  of 
oeean-vave  records  for  general  wave-height*  Diese  relations,  whioh  may  in 
part  be  deduced  from  oertain  assumptions  of  randomness  and  spectral  narrow- 
ness, appear  to  be  of  interest  both  as  explanations  of  previous ly-^nunn  wave- 
height  phenomena  and  as  possible  means  for  the  indirect  measurement  of  wave 
heights* 


2*  Harrow-spectrum  Theory  of  Have  Heights 

If  an  observed  ocean-wave  time  history  f(t),  measured  about  its  mean  value, 
be  regarded  as  a sample  from  a stationary  ergodio  random  process,  the  oon-^,% 
oept  of  an  associated  envelope  prooess  may  be  defined.  'Inis  has  been  done'  ' 
by  introducing  the  oonjugate  prooess  f*(t),  whioh  may  be  thought  of  as  being 
obtained  from  f (t)  by  shifting  each  of  the  angles  in  its  speotrum,  through  7T/2 
radians.  The  envelope,  defined  by  S(t)  = J ^ + [f*(t)3  will  be  a 

random  prooess,  the  ordinates  cn  the  graph  of  whioh  will  follow  some  proba- 
bility distribution  whose  nature  depends  upon  the  distribution  of  both  the 
ordinates  and  the  speotral  energy  of  the  original  prooess  f(t).  It  may  be 
noted  that  fi(t)  > f(t).  Pbr  the  limiting  o&ae  in  whioh  f(t)  = a sin  w t, 
we  have  f*(t)  & a eos  w t,  and  the  envelope  has  the  oonstant  value  given  by 
R(t)  ■ |4  Jbr  sufficiently  narrow-band  speotra,  E(t)  and  [R(t)3  ® may  be 
obtained  in  prinoiple  as  the  eleotronioally-filtered  outputs  of  ideal  linear 
or  square-law  deteotors* 

The  envelope  ourve  as  here  defined  may  be  shown  to  possess  the  property 
intuitively  expeoted  of  it  - namely,  that  it  pass  near  the  extreme  points 
on  the  original  ourve  - provided  the  relative  bandwidth  of  the  speotrum  of 
the  latter  is  small*  It  is  thus  not  unreasonable,  in  view  of  the  appe&ranoe 
of  ooean  swell  reoords,  to  expeot  the  distribution  of  the  wave  heights  on  the 
original  ourve  to  show  seme  agreement  with  the  distribution  of  the  ordinates 
of  itr  envelope* 


*k«e 


list  of  referenoes  at  end  of  report. 


3. 


Certain  as  swaptions  on  the  nature  of  the  ware-record  speotrum  have  been 
8 h wm ( - -- ) to  lead  to  a definite,  type  of  ordinate  distribution  for  the 

envelope  ourve  R(t)»  The  assvsnption  chosen  here,  which  is  in  part;  at 
least,  directly  verifiable,  is  that  expressed  by  the  oondition  that  f(t) 
be  a Qaussian  random  prooess (2X* (10), (ll).  Thig  assumption  has  the 
speoial  oonsequence  that  the  ordinates  of  the  original  ourve  have  a one- 
diasnsional  Gaussian  distribution^  That  is,  for  any  jr  and  for  any  dy  > o, 
tSai  urobaMll-by-  that  y~<  f tt ) < y + dv  is  given  by  l/(2  V 0~0^)Z  exp 
L~7*/(2<ro^)  J dy.  whhre  ; [_f(t)J  z,  the  bar  denoting  the  operation 
of  averaging*  A further  consequenoe  of  the  Gaussian  hypothesis  for  a 
stationary  ergodio  random  prooess  is  that,  in  a statistical  sense,  f(t).  and 
R(t)  are  desoribed  completely  by  the  unnormalized  correlation  function(l^) 

Y ( T ) * f(t)f(t  +i  )•  It  may  be  noted  that  Y (o)  s <r02,  a faot  whioh 
might  be  utilized  for  indirect  measurement  of  the  correlation  function. 

Under  these  assumptions  it  may  be  shown  that  the  ordinateB  of  the  envelope 
follow  a Rayleigh  distribution.  That  is,  for  any  r > o and  any  dr  > o, 
the  probability  that  r < r + dr  is  just  (l/ d02)rsxp  [ ”r2/^i  <r02)~j  dr, 

while  the  probability  that  R(t)  > r is  given  by 


rap  [-rV(?^o2)J 


0) 


It  will  be  seen  that  both  Rayleigh  and  Gaussian  distributions  depend  on 
the  single  parameter  <X02.  The  positive  quantity  o~Qt  the  standard  devi- 
ation, is  a measure  of  the  dispersion  of  the  ordinates  to  the  wave-rreoord 
olive  taken  about  their  mean  value  (assumed  here  to  be  zero).  Riysioally, 
for  purposes  of  measurement  it  may  be  interpreted  as  the  rco t-mean-s quare 
value  of  the  time-history  function  measured,  again,  from  its  mean.  Its 
square,  CTq2,  may  be  interpreted  as  the  average  total  power  (energy  per 
unit  time)  for  the  wave  speotrwa  of  long-ores  ted  linear  wave  theory. 


The  parameter  of  these  distributions  is  direotly  related  to  the  mean  value 
and  dispersion  of  the  envelope  distribution.  It  may  be  shown  that  the 
latter  parameters  are  respectively 


- <T0  (7T/8)*  = 1*253  C T0,  (2) 

and 

<3*s  r^o  r<4-TT)/2]i  = 0.685  cr0,  ...  (3) 

their  ratio,  oalled  the  relative  dispersion,  being 

[(t  ""Wf  = 0.523.  .. . (4) 

Sbrther,  a measure  related  to  distribution  asymmetry,  or  skewness,  previ- 
ously applied (I®)  to  wave-height  distributions,  is  found  to  have  the  value 


4. 


^ « +0*65  for  the  Rayleigh  distribution.  Finally,  the  mean,  R]^,  of 
the  highest  one-third  of  the  envelope  ordinates  is  given  by 

&1/5  9 c~0  ^(l*1  + 3 £ (2TT)~  - ^ «cp  (-t2/?.)  dt  j J 

...  (5) 

« 2.002  <r0  = yUR  4 1.143<TR  a 1.597 /*  R. 


Fbr  am  illustration  0?  the  way  in  whioh  the  Rayleigh  distribution  assigns 
different  probabilities  to  different  vertioal  intervals  on  the  wave-record 
chart,  one  say  rsfer  to  Figure  1,  where  the  straight  lines  represent  Rayleigh 
distributions.  The  horizontal  scale  has  been  laid  off  in  aooordanoe  with  (l), 
so  that  on  a straight  line  through  the  origin  a point  with  ordinate  r has  as 
abscissa  the  probability  that  R(t)&  r,  i.e.  the  peroent  of  the  time  Tike  env- 
elope ourve  is  no  higher  than  r.  The  mean  envelope  height  ,MR  is  then  the 
intercept  out  off  by  the  vertioal  line  erected  at  the  54.4#  level  (denoted 
by  the  oaption  "Mean"  on  Figure  1.)  The  standard  deviation  of  the  envelope 
is  given  by  the  difference  between  the  intercepts  of  the  line  in  question 
by  the  two  lines  labelled  "Mean  4 Dispersion"  (erected  at  the  level) 

and  "Mean".  The  standard  deviation  of  the  ordinates  on  the  original  wave 
reoord  having  the  given  envelope  distribution  would  be  cro  = /*r(2/tt)®* 

The  lines  on  the  figure  correspond  to  various  given  values  of  cr_*  the 
unite  being  the  same  as  those  for  the  envolope  ordinate  R(t).  It  is  seen, 
for  instance,  that  when  0~o  - 2.0  unite  we  have  95.6#  of  the  envelope 
ordinates  below  10.0  units  and  98.9#  below  12.0  units,  so  that  3.3#  of 
the  envelope  ordinates  are  between  10.0  and  12.0  units. 

The  interoept  out  off  by  the  line  ereoted  at  86.5#,  labelled  "R^ a"  in 
Figure  1,  is  equal  to  the  mean  of  the  highest  one-third  of  the  envelops 
ordinates. 


3.  Rmpirioal  Results  on  Estimated  Surfaoe  Wavs  Heists 

Previous  work^2)»^3^  has  led  to  results  on  the  distributions  of  wave 
heights  for  ocean  swell  in  shoaling  water.  These  heights  were  measured 
Aram  trough  to  orest  on  twenty-five  twenty-minute  subsurface  pressure 
records  and  then  modified  in  an  attempt  to  estimate  surface  wave  heists, 
i general  description  of  these  wave  rouorus  appears  ia  xaoxe  *.  uiq  re- 
sults obtained  whioh  are  of  present  interest  are  reoounted  below. 

The  dashed  ourve  in  Figure  2 represents  the  mathematical  distribution  type 
earlier  found  to  provide  a good  fit  to  the  estimated  ocean  surfaoe  wave- 
height  distributions.  This  model  distribution  is  the  member  of  the  family 
of  Gamma-type  distributions  with  the  skewness  coefficient  v 40.8.  It 
ts ay  be  ocmpared  with  the  dashed  straight  line  corresponding  to  the  Rijfleigh 
distribution.  There  has  also  been  drawn  on  this  figure  a full-line  ourve 
representing  the  Gsanaa-type  distribution  having  c<g  - 40.6,  which  is  the 
median  valua  fei*  the  twenty-one  linear  wave  records  examined  in  the  present 
study,  and  the  typical  value  found  earlier >3)  for  wind- generated  waves. 

All  three  distributions  ware  ohosen  to  have  the  same  mean  and  standard 
deviation,  so  that  their  shapes,  or  types,  oould  be  oorapared.  It  is  seen 
that  the  urevious  data  shew  little  avidanaa  for  th*  bypath1?*!” 

that  estimated  surfaoe  swell  heights  follow  a Rayleigh-type  distribution. 
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The  same  earlier  investigt/cion  established  the  dependence  of  the  wave-hei ght 
distribution  on  essentially  just  one  parameter.  This  scold  be  taken  ae 
either  the  mean  or  the  dispersion  of  the  wave-height  distribution,  sines  it 
was  found  that  these  two  parameters  were  quite  regularly  in  the  average 
ratio  = +0*53  - in  close  agreement  with  the  relative  dispersion  of 

the  Rayleigh  distribution* 

This  agreement  between  estimated  wave  heights  and  envelope,  both  in  distri- 
bution shape  and  parameter  ratios,  may  be  interpreted  as  evidence  that  the 
wave  records  studied  were  characterised  by  relatively  narrow-band  spectra. 
Ehe  regularities  in  wave-height  distributions  observed  earlier  were  used  to 
explain  the  regular  values  of  certain  nearly-oonstant  ratios  between  various 
wave-height  statistics  observed  by  other  investigators (17), (19), (20), (£2). 

It  is  now  soon  that  a more  fundamental  explanation  of  both  types  of  regu- 
larity is  provided  by  the  narrow-band  envelope  theory  oombined  with  the 
Gaussian  hypothesis  for  the  wave-reoord  ordinates* 


4*  Pressure  Wave  Ordinate  Results 

More  dirsot  verification  of  the  theory  presented  above  has  been  possible* 
Information  about  the  distribution  of  ordinates  on  the  pressure  wave-reoord 
ourve  was  obtained  in  January  1952.  This  distribution  was  found  to  be 
olosely  related  to  the  distribution  of  pr assure  wave  heights. 

Fbr  this  investigation,  the  Yave -Or dina ta  Distribution  Ivaalyzer  proposed  by 
and  cons trusted  under  the  supervision  of  P.B.  Snodgrass  was  used.  This 
electronic  equipment,  described  elsewhere (®7  in  detail,  receives  as  input 
the  varying  voltage  output  of  a pressure  meter,  or  of  a tape  or  other  record- 
ing of  this  voltage  si&ial*  Fbr  the  analysis  desoribed  here,  the  twenty- 
minute  pressure  reoords  previously (12), (13)  &£alyzod  a wave-wies  manner 
were  retraoed  by  a human  curve-follower,  producing  a replioa  of  the  original 
varying-voltage  signal,  oontaining,  however,  a relatively  small  amount  of 
added  "noise"  due  to  errors  in  tracing.  The  analyzer  performed,  in  effeot, 
a sampling  of  the  ordinates  to  the  voltage  ourve  every  0.1  second,  and  by 
means  of  ten  reoording  oounters  set  at  arbitrarily-chosen  voltage  levels, 
accumulated  the  total  number  of  sampled  ordinates  below  each  of  the  given 
levels. 

In  praotioe,  the  counter  voltage  levels  are  spaced  at  nearly  equal  intervals 
covering  the  entire  range  of  ordinates,  so  that  when  each  counter  reading 
is  divided  bytt>atotal  number  (approximately  12,000)  of  ordinates  sampleo, 
a definitive  set  of  points  on  the  cumulative  distribution  ourve  for  the 
ordinates  is  obtained* 

The  empirioal  ordinate  distribution  functions  for  the  twenty-one  reoords 
made  by  the  Mark  III  Pressure  Meter  (7)  are  plotted  in  Figure  3.  In  eaoh 
oaas,  the  ordinate  is  an  arbitrary  differential  pressure  level  measured 
in  ohart  divisions  and  the  abscissa  is  the  pereent  of  the  sampled  pressure- 
reoord  ordinates  whioh  do  not  oxoeed  that  level*  It  will  bo  seen  that  the 
plotted  points  appear  in  eaoh  ease  to  lie  on  an  approximately  straight  line* 
Sinoe  the  absoissa  soale  in  the  figure  has  been  taken  to  correspond  to  a 
Gaussian  probability  distribution  (s),  the  observed  results  indicate  a 


ORDINATE  DISTRIBUTION  FUNCTION 
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Gaussian  distribution  of  pressure-record  ordinates.  Similar  results  for 
wave  raoords  have  been  reported  by  J.E.  Meade  using  the  distribution 
density  analyzer  at  the  Naval  Hesearoh  Laboratory,  and  in  the  litera- 
ture . / » X IS J 


Bren  with  the  Gaussian  random  process  model  desoribed  above,  the  plotted 
points  in  figure  3 would  not  be  expected  to  lie  on  straight  lines  as  long 
as  the  ware  raoords  were  of  finite  length.  A measure  of  the  aotuai  sampl- 
ing fluctuation  to  be  expected  for  the  absoissa  of  any  one  point  on  these 
plots  may  be  oomputed  from  the  wave  speotrum  or,  equivalently,  from  its 
correlation  function.  The  amount  of  such  fluctuation  also  depends  upon 
the  temporal  spacing  between  oonsesutive  ordinates  sampled  from  the  wave 
reoord.  Vien  the  correlation  funotion  (or  the  speotrum)  for  the  wave  rec- 
ord is  also  known,  it  is  possible  to  assess  the  likelihood  that  the  sampled 
ordinates  come  from  a Gaussian  distribution  of  ordinates.  It  may  be  noted 
that  the  familiar  chi-square  test  often  used  for  this  purpose  is  based  upon 
independent  observed  random  variables,  i.e.  in  the  present  instance  upon 
ordinates  from  a wave  record  possessing  a completely  flat  and  broad  speotrum* 
Suoh  a situation  would  not  necessarily  be  expected  to  hold,  e.g.  for  a sampl 
ing  of  200-500  ordinates  from  a 20-30-minute  pressure  record  for  a typical 
correlation  function  indicating  a moderately  narrow  spectrum  oenterod  around 
10-15  seconds,  oven  if  the  sampling  is  done  in  some  random  manner*^-). 


The  result  of  the  analysis  of  Wave  Reoord  "B",  made  with  the  Mark  V Thermo- 
pile Sfavs  Meter , is  shown  in  Figure  4.  The  non-Gaussian  oharaoter  of 
the  ordinate  distribution  for  this  wave  reoord,  which  has  been  noted  by 
Birkhoff  and  Kotik'*',  is  shared  by  Wave  Records  "A*,  "C",  and  "D",  which 
are  the  other  three  Mark  7 records  among  the  original  twenty-five.  Owing 
to  the  non-linear  response  of  the  Mark  V instrument,  the  results  are  readily- 
interpreted  as  distortion  of  the  Gaussian  record  whioh  would  have  been  ob- 
tained with  the  linear  Mark  III  instrument.  Of  the  four  non-linear  time 
histories,  only  Wave  Reoord  "B"  was  analysed  by  the  Wave  Ordinate  Distribution 
Analyser.  For  the  other  three,  vicual  inspection  alone  was  sufficient  to 
show  considerable  positive  skewness  in  the  ordinate  distributions.  This 
oharaoteristio  was  also  noted  by  Dr.  Hiilip  Rudniek,  who  in  1951  analyzed 
Tfevs  Reoord  "B"  and  several  other  reoords  on  the  distribution  analyser  at 
the  TJ*S.  Naval  Eleotronios  Laboratory 


In  regard  to  the  slight  non-linearity  of  the  Mark  III  data,  no  interpreta- 
tion of  the  curvature  to  be  seen  in  some  of  the  plots  in  Figure  3 has  been 
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sis  ting  from  record  to  reoord,  would  be  oons  intent  with  an  expeoted  non- 
linearity^)* (1®/,  in  the  wave  reoords.  However,  on  the  reoords  analyzed' 
it  is  not  generally  known  whioh  is  the  direotion  toward  -the  crest  and  whioh 
is  that  toward  the  trough.  T£e  magnitude  of  the  depth,  varying  between  50 
and  86  feet,  at  whioh  the  pressure  reoords  'Were  made  is  probably  responsible 
for  there  being  no  disoernlble  relation  between  the  slight  degree  of  ourva- 
ture  in  the  plots  in  Figure  3 and  either  mean  wave  height  or  depth  of  record- 
ing. 


The  twenty-one  plote  in  Figure  3 not  only  provide  evidence  for  the  existenoe 
of  underlying  Gaussian  ordinate  distributions,  but  the  individual  plots  may 
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also  be  used  to  estimate  the  parameters  of  those  distributions.  Various 
general  methods  of  estimating  the  ordinate  dispersion  and  the  ordinate 
mean  /^0  were  employed,  in  order  to  investigate  both  the  validity  of  the 
envelope  theory  and  the  aocuraoy  obtainable  in  rapid  graphioal  estimation 
of  these  parameters.  Only  two  of  these  methods,  referred  to  as  the  methods 
of  "direot  computation"  and  "graphical  estimation",  are  discussed  her©  and 
in  Section  6.  Chief  interest  in  the  present  investigation  was  in  the  rela- 
tion between  the  ordinate  dispersion  and  the  distribution  of  wave  heights. 
For  many  purposes  on  estimate  of  the  mean  ordinate  A*0  may  be  required,  aB 
when  measurements  are  to  be  made  from  "still-water  level",  as  assumed  in 
the  theory. 

The  ordinate-distribution  parameters  for  each  wave  reoord  were  computed 
direotly  from  the  statistical  frequenoy  distribution  composed  of  finite- 
width  olasses  defined  by  the  counter  levels  and  the  eounter  reading. 

Thie  computation,  to  whioh  a correction  for  the  finite  average  class- 
width  was  applied,  is  regarded,  for' purposes  of  comparison  with  the  theory, 
as  yielding  the  iaost  reliable  estimates  of  <rQ  and  Q,  whioh  are  taken 
aa  standards  for  comparison  with  wave-height  statistics  and  other  subse- 
quent estimates  of  these  distribution  parameters.  The  direotly-computed 
values  of  <0"o  appear  in  Table  I.  Because  the  set  of  counter  levels  used 
in  analysing  the  data  rarely  included  the  entire  range  of  ordinates  present 
on  the  wave  reoord,  extension  of  the  olasses  provided  by  the  oounter  levels 
was  made.  This  was  dene,  so  as  to  give  ten  classes,  by  assuming  that  each 
ond-o las s required  to  encompass  one-hundred  peroent  of  the  ordinates  would 
have  the  same  width  as  the  olass  adjaoent  to  it. 

The  full  lines  appearing  on  Figure  3 represent  Gaussian  distributions, 
their  intercept  and  slope  being  in  each  oaBe  the  direotly-computed  mean 
and  standard  deviation.  The  fit  of  eaoh  of  the  lines  to  the  pressure 
time-history  ordinate  data  is  seen  to  be  reasonably  good. 


6.  Pressure  Ordinates  and  ffava  Heights 

Aa  & check  on  the  application  of  the  narrow-band  envelope  theory  to  pressure 
rooords  of  ooean  swell,  the  relation  between  the  observed  ordinate  dispersion 
and  mean  treugh-tc-orest  wave  height  was  examined.  For  this  purpose,  a wave 
height  on  the  twenty-one  pressure  records  was  oouated  and  measured  between 
eaoTflai n Imurn  (trough)  and  the  immediately  following  maximum  (orest).  The 
mean  value  of  these  maximua-minimiaa  pressure  differences  is  referred  to  as 
the  mean  wave  height  [Agu  The  computed  values  of  this  quantity  are  shewn  in 
Table  I.  The”  relation  between  A1 2 and  the  ocmputed  value  of  <rQ  is  shown  in 
Figure  5,  together  with  the  theoretioal  straight  line  for  the  mean  2 
of  the  double -envelope,  whose  equation,  obtained  by  doubling  both  members 
of  (2),  i's  2 A*  p s 2 0~0  £ /2)®»  Agreement  with  this  line  would  be  ex- 

pected only,  or  oourse,  if  measurements  were  made  on  the  true  envelops  of 
an  infinite- length  wave  reoord  satisfying  the  Gaussian  process  hypothesis. 

A straight  lias  fitted  to  the  data  is  seen  to  be  approximately  15jJ  below 
the  envelope  line.  The  slope  of  the  1 sas-c-squares  line  through  the  origin 


shown  in  Figure  5 is,  in  f&ct,  0.850  times  the  slope  of  the  envelope  line. 
This  systematic  deviation  of  the  data  is  in  the  expected  direction,  since 
even  the  extremes  on  the  wave-record  curve  cannot  lie  outside  the  envelope 
ourve.  The  magnitude  of  the  deviation  would  be  expected  to  vary  from  one 
wave  record  to  another  in  accordance  with  the  relative  bandwidth  of  the 
speotrum.  Present  indications  are  that  a rough  estimate  of  this  spectral 
characteristic  will  provide  enough  information  to  eliminate  a large  part 
of  the  average  discrepancy  between  the  points  and  the  line  in  Figure  5. 

3ased  on  the  present  data,  the  empirical  relation  between  me&r?  wave  height 
and  ordinate  dispersion  becomes 

H = 2(0.85)  3-0  ( Tc/2)&  ...  (6) 


Substituting  this  empirioal  relation  for  the  relation  (2),  we  may  replace 
the ' straight  lines  of  Figure  1 by  those  in  Figure  6,  whioh  shows  the  result- 
ing semi-empirical  set  of  wave-hei$it  distributions  for  various  ordinate 
dispersions.  Given  the  ordinate  dispersion  <TQt  the  corresponding  line  on 
Figure  6 represents  the  expeoted  distribution  of  wave-heights  given  by  the 
Rayleigh  distribution  whose  parameters  are  obtained  tram  the  line  fitted 
to  the  data  in  Figure  5. 

As  regards  the  so&tter  of  the  points,  the  significant  fact  wh,ioh  K,gure  5 
illustrates  is  that  the  sample-length  and  spectral  bandwidth  of  the  pressure- 
wavs  records  analysed  wers  of  suitable  magnitudes  to  oreate  a fairly  olose 
relationship  between  the  measured  estimates  of  mean  wave  height  and  ordin- 
ate dispersion.  In  other  words  & knowledge  of  the  measured  ordinate  r.m.s. 
value  (or  the  total  spectral  power)  for  these  wave-record  samples  is  essen- 
tially equivalent  to  a knowledge  of  the  measured  moan  wave  height,  within 
5%  over  two-thirds  of  the  time.  Fbr  longer  wave  reoords  it  may  be  expeoted 
that  the  peroent  of  soatter  in  Figure  5 will  be  reduoed.  As  has  been  seen, 
the  mean  wave  height  essentially  determines  the  distribution  of  all  wave 
heights.  Hiiat  the  present  study  shows  is  that  sampling  fluctuations  do  not 
seriously  alter  the  effective  validity  of  the  theoretical  relationship  when 
20-minute  pressure  reoords  sure  used. 

The  most  important  implication  of  this  feet  is  that  for  such  wave  reoords, 
n wave -by-wave  analysis  for  wave  heights,  a prooedure  relatively  difficult 
to  mechanize,  may  be  replaced  by  an  ordinate  analysis  for  O'..  Since  the 
ordinate  dispersion  O',  can  be  oomputed  as  an  r.m.s.  value,  the  latter  task 
may  be  aoocaaplished  in  principle  by  means  of  a direot-reading  eieotrioal 
meter  provided  with  a sufficiently  long  time -constant. 

Under  same  oiro  urns  tone  es  the  measurement  of  <T0  s { [f(t)3  2 may  be 
l«es  convenient  than  that  of,a0  s |f(t}( , the  mean  absolute  deviation  of 
the  ordinates  to  the  wave  reoord  (measured  Aron  the  mean  ordinate).  Eleo- 
trio&l  measurement  of  ao  is  possible  in  principle  by  means  cf  a rectifying 
circuit.  Sinoe  for  a Gaussian  distribution  we  have 

c r0  s a0  (-»r/8)*,  ***  (7) 

the  mean  envelope  height  is  given  by  the  theoretical  relation 
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*nd  the  mean  wave  height  by  the  empirioal  relation 

/*H  s 2(0.86)  a0  (77/2)  ...  (9) 

The  sampling  stability  of  the  ratio  t^/ CT-,  whose  expected  value  is 
(2/rr )a«  0.80,  given  by  (7),  is  indioatea  in  the  results  published  by 
SsiwellC^8)  for  wave  reoords  slightly  less  than  six  minutes  in  length. 

The  values  were  reported  to  be  within  the  limits  0.78  - 0.84,  a range  in 
which  a reduotion  could  be-  expeotod  if  twenty^  or  thirty-minute  wave  reo- 
ords were  taken. 

Obviously,  the  relations  developed  here  between  various  distribution 
parameters  mi  git  be  made  the  basis  of  a number  of  graphs  similar  to  Fig- 
ures 1 and  6,  showing  the  distribution  of  either  the  theoretioal  envelope 
R(t)  or  the  wave  heights  (IL  ) for  various  values  of  any  one  of  the  follow- 
ing parameters t Ordinate  dispersion  (O'),  ordinate  mean  absolute  deviation 
(ao)f  mean  wavs  height  (^g),  wave-height  dispersion  (ff"H),  mean  highest 
one-third  wave  height  envelope  mean  (A*®),  envelope  dispersion 

(<Tg)  or  envelope  mean  highest  one-third  (R^yw)*  Suoh  graphs  would  be  of 
varying  interest,  depending  upon  whioh  parameter  was  oonsidered  as  being 
the  basic  measured  one. 

Sinoe  the  most  fundamental  parameter  is  the  ordinate  dispersion  <TQ  (or  the 
total  epeotral  power  <0'(,.^),  we  summarise  in  Table  II  the  relations  between 
the  parameters  by  giving  the  ratio  of  eaoh  to  <T0.  The  ratio  of  any  other 
pair  of  parameters  may  be  found  ft*cm  these  ratios  by  performing  a single 
division.  The  ratios  in  Table  II  may  be  used  to  replace  the  parameters 
marked  oh  the  straight  lines  in  figures  1 and  6 by  any  other  parameters, 

••S*  /* H»  Hyfc»  or  *0* 


Table  II 


Ratios  of  Various  Distribution  Parameters  to  Ordinate  Dispersion 


parameter  | 

*0 

/*R 

^R 

ESS! 

^ H 

EiA 

I 

(77/2)* 

[(4-TT)/^ 

(of.pg.i) 

1,7 

ES3 

rfl/s) 

(parameter)/  ^ 

K 

■HH 

mOU 

Wo  / 

1 

miH 

1.253 

0.666 

2.002 

2.130 

HU 

3.403 

In  making  use  of  these  ratios  it  must  of  oourse  be  remembered  that  they  de- 
scribe theoretioal  probability  distributions,  and  oan  only  be  expeoted  to 
approximate  olossly  the  oonputed  statistics  of  an  empirioal  frequency 
distribution  if  the  latter  are  based  on  a sufficiently  long  sequenoe  of 
data,  i.e.,  a sufficiently  longwave  reoord.  The  magnitude  of  the  sampling 
fluctuations  in  certain  statistics  for  twenty-minute  wave  reoords  has  been 
indioated  by  results  in  previous  reports (12), (13), 
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6*  Simplified  Estimation  of  Pressure  Ordinate  Dispersion. 

The  near-reoti linearity  of  the  plots  in  Figure  3 makes  possible  for  eaoh 
wave  record  a convenient  graphical  estimate*  of  ths  mean  ordinate  f*  and 
the  root-mean-square  (or  standard)  deviation  'j'q  fro a the  mean  ordinate. 

The  first  is  read  directly  off  the  graph  as  tho  intercept  at  the  509?* 
point  on  a visually- fitted  line.  The  second  is  "taken  as  one-half  the 
difference  between  the  intercepts  at  the  S4. and  15.9#  points,  a measure 
of  the  slope  of  the  line. 

The  degree  of  association  between  the  graphical  estimate  of  0”o  and  the 
mean  wave  height  is  shewn  in  Figure  7.  The  graphical  estimate  of  <TQ  ob- 
tained from  the  visually- fitted  line  on  Gaussian  probability  paper  appears 
to  be  equivalent,  if  multiplied  by  a constant,  to  the  mean  pressureHrave 
amplitude  within  about  12#  in  nearly  all  oases.  In  fact  this  graphical 
method,  carried  out  before  the  dir eot-o amputation  method  was  employed, 
appears  to  be  almost  as  good  as  the  latter,  as  can  be  seen  by  a comparison 
of  Figures-  S and  7. 

It  is  evident  that  for  wave  systems  having  ordinate  distributions  suffi- 
ciently close  to  the  Gaussian  type,  the  oomplete  sets  of  points  shewn  on 
Figure  3 should  not  be  neoessary  for  obtaining  reasonably  good  estimates 
of  the  slopes  crQ.  Far  example,  it  might  be  possible  to  use  two  oounters 
set  at  ordinate  levels  yQ  and  yp  ohosen  in  aooordanoe  with  the  expeoted 
ordinate  dispersion,  i.e.  the  average  wave  height  to  be  encountered. 

A two-counter  analyser  might  be  suitable  for  field  installations  in  whioh 
eifcher  a stamped  or  a photographio  record  of  counter  readings  made  at 
regular  intervals  would  be  obtained  as  the  only  data. 

To  estimate  the  ordinate  dispersion  from  two  oounter  readings,  one  would, 
as  in  Figure  3,  plot  the  two  oounter  levels  y0,  yi  against  the  two  relative 
oounter  readings,  P(y0),  P(yi),  obtained  by  dividing  the  actual  oounter 
readings  by  the  total  count  Twhich  might  conveniently  be  made  12,500  or 
20,000  0.1-seoond  time  units).  These  two  points  determine  a unique  line 
on  whioh  a oertain  increment  measured  on  the  vertioal  soale  is  out  off  by 
ttxe  6005  and  the  84.  ]£  points:  The  quantity  <r^’,  the  desired  estimate  of 

"is  is  this  vsrtioal  inoreaent,  whioh  may  also  be  oomputed  numerically 
by  expressing  the  slope  of  the  line  analytical?^’  as  a function  of  yQ,  yj, 
P(y0),  and  P(y^),  and  making  use  of  a tf-.ble  of  the  Gaussian  probability 
function  or  the  error  function* 

The  latter  pro cs dure  was  adopted  for  computing  the  estimates  <T0',  whioh 
are  plotted  versus  the  direotly-ocmputed  C0  in  Figure  8.  To  obtain  these 
estimates  for  eaoh  wave  reoord,  y0  and  y^  were  seleotod,  as  described  below, 
from  the  ten  oounter  levels  used  in  analysing  the  data.  It  will  be  seen 
that  the  estimate  of  ordinate  dispersion  obtained  from  the  two- counter 
data  is  relatively  olpse  to  that  obtained  using  all  ten  oounters.  Only 
ones  among  these  wave  reoorda  does  it  differ  by  more  than  ten  peroent, 
and  over  four-fifths  of  the  time  the  relative  difference  was  less  than 
five  peroent* 

The  estimation  procedure  may  be  further  simplified  by  eliminating  one  of 
Ihe  two  counters  if  y0  oan  be  taken  as  the  mean  ordinate  or  the  ordinate 
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corresponding  to  the  sero  differential-pressure  level,  whose  location  may 
ba  pre-detemined  for  a given  instrument  installation.  If  this  is  don*, 
then  P(y0)  will  be  0*50,  and  one  of  the  two  plotted  points  is  fixed.  Bis 
location  of  the  other  point  (y^,  P(y^)  ) determines  the  estimate, 
which  will  be  obtained* 

figure  9 shows,  for  various  oounter  levels  y-,s  the  carves  giving  the  estim- 
ated value  <r0 " as  a funotion  of  the  single  counter  reading  F(y^)  expressed 
in  peroent  of  the  total  oount.  For  example  if  the  oounter  is  set  at  4 feet 
and  reads  £4.3$,  the  estimate  <T0 B is  2*5  feet,  fbom  this  value  it  may  also 
be  Inferred  that,  e.g. , HWj  is  2(0*86)  (4.0)  - 6*8  feet,  etc.  As  explained 
earlier,  use  may  be  made  or  the  ratio  between  0"o  nnd  any  other  parameter 
of  interest  to  make  the  vortical  aoale  in  figure  S dlraot-reading  for  that 
parameter  - «*g*  /Wg,  H^,  eto. 

While  the  average  values  of  the  estimates  <T0"  are  indicated  by  Figure  9, 
the  statistioal  fluctuation  of  the  estimates  based  on  the  one-counter 
analysis  of  a twenty -minute  reoord  may  be  visualised  from  Figure  3.  Bis 
two  dashed  lines  on  the  plots  in  Figure  3 have  been  drawn  to  correspond 
to  the  upper  and  lower  10$  limits  for  the  error  in  predicting  cr_  from. 
individual  plotted  points  when  the  mean  ordinate  M ~ is  known*  The  estim- 
ate of  <Te  obtained  in  this  way  frcm  any  one  point  Is  within  410$  of  the 
value  computed  direotly  from  all  ten  points  if  and  only  if  that  point  lies 
in  the  two  narrow  seotors  between  the  taro  dotted  lines.  It  can  be  seen 
that  percentage  errors  are  apt  to  be  high  for  single  points  vhioh  are 
either  too  near  the  center  of  the  distribution  or  else  too  far  out  on  its 
extremes*  There  are  indications  that  the  most  reliable  estimates  of  CT0 
are  obtained  when  the  ordinates  y^  and  y0  are  ohoaen  symmetrically  about 
the  mean  at  s.  distance  of  about  1*6  0~0,  Thors  positions  in  the  daussian 
distribution  of  crdinates  correspond  to  oounter  settings  at  the  93/5  and 
the  7$  points,  respectively*  The  dotted  vertioal  line  in  Figure  8 allows 
this  setting  of  y-j,  to  ba  read  off  direotly  for  any  given  ordinate  dispersion 
<T0«  The  oounter  readings  used  for  the  estimates  <TC»  shown  in  Figure  8 
were  chosen  to  be  those  most  nearly  at  thssa  vertioal  levels* 

Although  CT0  is  initially  unknown,  a relativsly  small  amount  of  preliminary 
observation  would  allow  the  oounter  to  be  set  for  further  observation  near 
the  optimum  level  +1*6  crot  or  the  96$  point,  in  terms  of  the  oounter  read- 
ing* for  a permanent  installation  yi  ocuid  be  automatically  adjusted  #t 
intervals,  if  naoassary*  Fbr  example,  the  oounter  level  may  be  shifted 
every  half-hour  to  41*6  times  the  previous  estimate  of  c T0,  thus  keeping 
approximately  in  step  with  the  clew  changes  in  the  total  wave-speotrum 
power*  If  such  shifts  in  y^  were  made  from  time  to  time  to  take  aooount 
of  changing  wave  oenditiens,  it  would  be  necessary.  of  oourse,  to  reoord 
each  new  level  y^» 

In  the  absence  of  self-adjusting  counter  levels,  one  or  two  fixed  levels 
oould  be  chosen  as  often  as  feasible  in  aooordanoe  with  the  expected  values 
of  <T0  for  the  tiv.ve  conditions  of  interest  most  likely  to  ooour  over  an 
observation  interval  of,  perhaps,  days  or  weeks*  In  this  oase  a knowledge 
of  the  wave  climate  for  the  plaoe  and  season  of  installation  oould  be  used 
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as  & basis  for  selecting  the  fixed  ordin&te  levels.  This  type  of  wave- 
olioate  information  should  be  available  after  sufficient  data  taken  in  this 
nay  at  the  installation  site  has  accumulated.  In  making  use  of  suoh  informa- 
tion a oompromlse  would  be  found  necessary  between  (l)  the  loss  of  aoour&oy 
for  high  wave  reoords  when  the  oounter  level  is  oho sen  too  low,  and  (2)  the 
loss  of  nearly  all  information  from  low  wave  reoords  when  the  oounter  level 
is  ohosen  too  hig^i. 

As  an  example  of  the  application  of  the  fixed- level  procedure  to  a given 
wave-meter  installation,  the  following  hypothetical  situation  may  be 
imagined.  Suppose  that  the  twenty-one  wave  reoords  B,  F,  ...,  T had  been 
obtained  at  a single  location  (instead  of  four)  over  a single  season  of 
the  year  (instead  of  all  seasons).  If  this  data  were  oonsidered  as  histori- 
cal wave-olimate  data  for  suoh  an  installation,  the  anticipated  average  valuo 
of  cro  would  be  approximately  6.0  ohart  divisions  above  the  mean.  The  low 
wave  reoords  show,  however,  that  unless  the  oounter  level  is  made  as  low  as 
about  4.0  units  above  the  mean,  no  information,  exoept  an  upper  limit  to  the 
ordinate  distribution,  would  be  obtained  in  some  oases. 

A rough,  probably  conservative,  test  of  this  one-oounter,  fixed-level  method 
of  estimating  CTQ  was  oonduoted  by  o hooting  the  level  y^  to  be  the  level  of 
the  oounter  which,  in  the  original  analysis  of  the  present  data,  was  set 
nearest  to  4.0  ohart  divisions  above  the  mean.  As  a result,  for  either 
(l)  low  waves  suoh  as  in  Wave  Heoord  "0"  or  (2)  high  waves  suoh  as  in  Wave 
Reoord  *V",  the  point  ohosen  is  f hr  from  the  93#  level,  and  the  estimated 
slope,  cr  ",  is  subjeot  to  relatively  large  errors.  These  will  result  from 
either  (1)  division  by  a small  number  or  (2)  multiplication  by  & relatively 
instable  observed  value  far  out  in  the  tail  of  the  ordinate  distribution. 

From  the  oounter  levels  shown  in  Figure  3 it  may  be  seen  that  most  of  the 
points,  the  distance  of  whose  ordinate  above  the  mean  is  nearest  4.0  ohart 
divisions,  would  fall  within  a tlOjC  error  band.  The  values  of  the  estimate, 
crQmt  were  oomputed  analytically  and  plotted  in  Figure  10  against  the  values 
of  <ro  already  oomputed  by  using  the  data  from  all  ten  oomters.  The  two 
outer  lines  in  the  figure  indicate  the  8 % error  limits. 

The  advantage  of  an  approximate  knowledge  of  the  93jt  ordin&te  level  is 
shown  by  & comparison  of  the  estimates  plotted  in  Figure  10  with  those 
plotted  in  Figure  11.  The  latter  estimates  ” were  based  on  a knowl- 
edge of  the  mean-ordinate  level  and  were  obtained  from  the  reading  of  the 
oounter  nearest  the  $$%  level  (i.e.  one  of  the  oounters  used  in  ocmputing 
the  two-oounter  estimate  of  the  ordinate  dispersion). 

Figures  8,  10  and  11  indioate  that  satisfactory  estimation  of  the  r.sua. 
value  (and  henoe  mean  wave  height,  eto.)  of  a wave  reoord  may  be  aohleved 
by  the  use  of  only  one  or  two  ordinate  oomters.  Suoh  a routine  estimation 
should  provide  useful  information  about  wave  reoords,  if  used  over  & suit- 
able length  of  time  in  a nisnber  of  installations.  If  a sequence  of 
estlmatss  <TQ"  were  made,  say,  at  ons-half  hois*  intervals,  the  results 
would  bs  to  some  extent  se If -ohe eking,  in  that  values  dlffereing  markedly 
from  the  general  trend  oould  be  interpreted  as  having  bean  overly  influenced 
by  random  factors.  Thus  significant  trends  in  <ro  reflecting  100-200JC 
ohanges  in  wave  energy  would  not  be  readily  obsoured  by  an  occasional 
random  lOjf  change  observed  in  CT0". 
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Oo (Min -ware  reoord  ordinates  are  found  to  have  approximately  the  follow- 
ing properties  of  a narrow-speotrum  Gaussian  random  prooec ■ i (l)  Pressure 
wave-reoord  ordinates  follow  closely  the  Gaussian  distribution*  (2)  Estim- 
ated surfaoe  wave  heights  follow  olosely  the  Rayleigh  distribution.  (S)  Est- 
imated surfaoe  ware -height  distribution  parameters  are  related  to  one  another 
in  the  same  way  as  are  those  of  the  Rayleigh  distribution,  all  being  olosely 
proportional  to  the  pressure  ordinate  dispersion  C~Q,  ($•)  Pressure  wave- 
height  mean  and  pressure  wave-ordinate  dispersion  are  olosely  proportional* 

fbr  pressure  wave  reoords,  wave-height  mean  is  found  to  be  approximately 
1.70/  times  wave-ordinate  dispersion.  Graphloal  estimation  of  wave- 

ordinate  dispersion  is  found  to  be  a reasonably  aeourate  method  of  estimat- 
ing wave-height  mean  for  twenty-minute  pressure  reoords  of  ocean  swell. 

A fixed, suitably-situated, single  point  on  the  ordinate  distribution  ourve 
is  fo^md  to  yield  a reasonably  aoourate  estimate  of  ordinate  dispersion  if 
the  "still-water  level"  is  known,  while  the  use  of  one  or  two  variable 
points  gives  improved  results. 

A systematic  program  of  ocean-wave  measurement  utilising  these  simplified 
forms  of  the  Wave-Ordinate  Distribution  Analyser  may  be  expeoted  to  provide 
useful  continuous  wave-height  information  in  immediately -available  form. 
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